MINIMALITY VIA SECOND VARIATION FOR MICROPHASE SEPARATION OF 

DIBLOCK COPOLYMER MELTS 



VESA JULIN AND GIOVANNI PISANTE 

Abstract. We consider a non local isoperimctric problem arising as the sharp interface limit of the 
Ohta-Kawasaki free energy introduced to model microphase separation of diblock copolymers. We 
^~»^ i perform a second order variational analysis that allows us to provide a quantitative second order min- 

^s^J ■ imality condition. We show that critical configurations with positive second variation are indeed strict 

£H ■ local minimizers of the nonlocal perimeter. Moreover we provide, via a suitable quantitative inequal- 

ity of isoperimetric type, an estimate of the deviation from minimality for configurations close to the 
minimum in the L topology . 

o 

m 

1. Introduction 

< 

In this note we are interested in performing a second order analysis for a nonlocal isoperimetric 
problem arising as a variational limit of the Ohta-Kawasaki functional introduced for a density functional 
theory for microphase separation of A/B diblock copolymers. 

Among the several mean field approximation theories proposed to model the phase separation in 
diblock copolymer melts, the one derived by Ohta and Kawasaki in [15] turns out to be one of the most 
promising from the mathematical point of view. Let il C M. n be a domain representing the volume 
occupied by the polymeric material. The free energy can be written as a nonlocal functional of Cahn- 
Hilliard type as 

(1.1) £ s {v)=e I \X/u\ 2 dx + - I {u 2 - l) 2 dx + 70 / / G(x,y)(u(x) - m)(u(y) - m) dxdy, 
■ Jn £ Jn Jn Jn 

where u G H (O) represents the density distribution of the monomers forming the copolymers, G is the 



Green's function of the Laplace operator with Neumann boundary conditions and 

u dx 



n 

is the difference of the phases' volume fractions. The formulation (jl.l[) was first introduced in [14] and 
for a derivation of the Ohta-Kawasaki density function theory from the self-consistent mean field theory 
we refer the reader to [6] and the references therein. 

As pointed out in [B] , depending on the molecular structure of the polymers there are several regimes 
of phase mixture. Nevertheless the presence of an observable phase separation occurs in the so called 
intermediate or strong segregation regimes, wherein the domain size is much larger than the interfacial 
length. This suggests that most of the features of the model can be described, from a mathematical point 
of view, by looking at the sharp interface limit of £ £ as the thickness e of the diffuse interface tends to 
zero. Thus we are lead to study the minimizers of the following energy functional, that arises as the 
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r-limit of £ £ in the L 1 topology, 

(1.2) £(u) = ^\Dum)+ 7o f f G(x,y)(u(x)-m)(u(y)-m)dxdy 

6 Jn Jn 

where u : — > { — 1, 1} is a function of bounded variation and |Du|(f2) denotes its total variation in tt. 
The functional (jl.2p can be described in a more geometric way that turns out to be more suitable for our 
analysis. Indeed identifying the function u with the set E = {x G il \ u(x) = 1} and using the properties 
of G(x,y), we can rewrite £ (it) as 

(1.3) J{E) = P n (E) + 7 / \Vv E \ 2 dx, 

Jn 

where Pq(E) stands for the perimeter of E in Q and ve is the solution of 



(1.4) 

with the positions 



-Ave = ue — m in 
^ = on dCl, 



ue dx. 



ue =XE~ Xn\E and m = I 

Jn 

We are thus lead to study the variational problem 

(1.5) min | E C O, J u E dx = 

We explicitly observe that the condition J^UEdx = m is nothing but a volume constraint since \E\ = 
i(m+|0|). 

The problem (jl.5p has been presented in [7] as a mathematical paradigm for the phenomenon of 
energy-driven pattern formation associated with competing short and long-range interactions. Recently 
an increasing interest has been devoted to the second order variational analysis of energy functionals 
which exhibit this competing behaviour. Indeed it has been successfully applied to prove stability and 
minimality criteria in several contexts (see for instance [HE]). Directly related to our problem is the work 
by Ren and Wei and by Choksi and Sternberg (cfr. [TCI [T71 [23 HH [201 In a series of papers, as a first 
step toward the validation of the conjectured periodicity of global minimizers of (|1.3[) , they calculate the 
second variation of (|1.3[) at a critical configuration and, among several other applications, they construct 
examples of periodical critical configurations and find conditions under which their second variation is 
positive definite. 

A different approach is taken in [1], where the authors consider the problem of minimizing the 
functional (|1.3|) in the periodic case, i.e., when Q is the n-dimensional fiat torus. They prove the local 
minimality of critical configurations which second variation is positive definite. In this direction goes also 
our present work, whose aim is to prove that a similar minimality criterion holds true for the problem 
(|1.5[) . Indeed we prove in Theorem 12.31 that if E is a sufficiently regular critical point of (|1.3p such that 
the quadratic form d 2 J(E) associated to the second variation of J at E is positively defined, then there 
exists a constant Co > such that 

J(F) > J(E)+c„\FAE\ 2 

for any admissible set F C f2 sufficienlty close to E in the L 1 topology. This in particular implies that E 
is a local minimizer of J and in addition provides a quantitative estimate of the deviation from minimality 
for sets near to E. We remark that our result holds also in the case 7 = and thus we cover the case of 
volume constrained isoperimetric problem in a regular domain Q. 
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Some comments are in order on the differences between the periodic case (studied in [T]) and the 
Neumann boundary case. Indeed if on one side working with the Neumann setting dispenses us from 
several technicalities needed to deal with the translation invariance of the functional J, on the other side 
new delicate arguments need to be introduced to deal with problems which arise when E touches dVL. We 
remark that in [1] the Neumann boundary case was considered only under a rather restrictive assumption 
that the set E does not intersect the boundary of SI. 

Finally we outline the structure of the paper. In the next section we introduce the notation and 
present the framework in order to precisely state our main results. In Section|3]wc discuss some regularity 
results for oj-minimizcrs of the area functional and in particular a stability result for the regularity (cfr. 
Theorem I3.5[) which will play an important role in the proof of the main theorem. Section 4 is devoted 
to the lengthy calculations of the second variation formula for regular sets satisfying the orthogonality 
condition on <9S1. The proof of the main result, Theorem 12. 3[ is divided in sections 5 and 6. The scheme 
follows a well established path (see for instance [HIS]). First we use the general second variation formula 
from Proposition 14. II to prove the local minimality among regular sets which are close to the critical one 
in the VF 2 ' p -topology. The final result is then proved by contradiction using a penalization argument and 
exploiting the regularity theory of w-minimizers. 



In this section we set up the basic notation, recall some preliminary results and present the statements 
of the main results. Throughout the paper S7 C K n is assumed to be a bounded domain with C 3,a - 
boundary, for some a > 0. 

We say that a set E C SI is C k ' a regular, where k < 3, if its relative boundary M 



is a C fc,a -manifold with or without boundary if dE n 9S1 ^ or dE n <9S1 = respectively. Suppose E is 
C 1 regular and let X be C 1 -vector field in SI which satisfies 



2. Preliminaries and statement of the results 



(2.1) 



M := dE n SI 



(2.2) 



X(x)eT x (dQ) for every x G dSl, 



and 



(2.3) 




where vm is unit normal of M. We define the associated flow $ : Q, x (— to, to) — s- SI by 




$(x,0) 



We then define the first variation of (|1.3[) 



at E with respect to the field X (or flow $) by 




t=o 



and the second variation by 




t=o' 



where E t = <$>{E 7 t). 

Formula for the first variation of (|1.3[) is well known. 
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Proposition 2.1. Let E <Z Q, be a C 1 -regular set with ue dx = m. Suppose that X satisfies (|2.2|) and 
()2.3|) . The first variation of (|1.3j) at E with respect to X can be written as 

dJ(E t 



t' I _ / l TJ - _ i A~,„. _\/v 1 I / ,,*\ ^oj«-2 



#M + 4 7 ^)(X, l/ A/ > cM"" 1 + / (X, I/*) d?T 

m JMnon 



Here v* is the outward unit co-normal of M n 90, i.e., normal to M PI 90 and tangent to M , and uq is 
the unit normal o/O. 

We say that C 1 regular set E C O is critical if the first variation is zero for every C 1 -vector field X 
which satisfies ([2~2]l and ((2~3")) . 

If we choose the vector field X such that it is compactly inside O, spt(X) (s O, then the last term in 
the first variation is trivially zero. Therefore if E is critical it satisfies the Euler-Lagrange equation in a 
weak sense 

(2.4) H M + 4<yv E = A on M, 

where the constant A is Lagrange multiplier associated to the volume constraint. We may further deduce 
that also the second term in the first variation vanishes 

(X, v*) dU n - 2 = 0. 

Mndn 

This and (|2.2|) imply that M is orthogonal to 90 on M n 90. These observations motivate us to define 
critical sets for (11.311. 



Definition 2.2. We say that C 1 regular set E C O is a regular critical set if it satisfies (|2.4[) in a weak 
sense and meets 90 orthogonally, i.e., 

(2.5) (u M , m) = on M n 90. 

We may use standard regularity theory for elliptic equations, as in [T], to deduce that every regular 
critical set is a C 1,a -manifold with boundary. Since ve solves (|1.4|) we obtain that ve G C 1: ^(0) for every 
/? G (0, 1). Standard Schaudcr estimate then implies that M is in fact C 3,,3 -mamfold with boundary for 
every f3 <G (0, 1). 

It is well known that the second variation of (|1.3[) at a regular critical point has a special structure 
and it can be written as a quadratic form. To present the formula we remark that the non-local part of 
(|1.3jl can be written more explicitly in terms of E by using Green's function. Let G : O x O —> M. be the 
Green's function with Neumann boundary conditions 

-A y G(x, y) =S X - p- in O, 

= on 90, 
j n G(x,y)dy = 0. 

The functional (| 1 . 3|) can be written as 

J(E) := P Q (E) + 7 f [ G(x, y) u E (x) u E (y) dydx. 

Derivation of the second variation formula for critical sets can be found in [7| and [35]: as it was 
noted after Remark 2.8 in [J. To present the formula let E c O be a regular critical set and denote its 
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relative boundary by M. The second variation of (|1.3|) at E with respect to X can be written as 

d 2 J{E)[<p] = f (\D T ^\ 2 - |5 M |V 2 ) dH^ 1 - f B m (y M ,vu)v 2 dU n - 2 
Jm JMndn 

(2.6) +87/ [ G(x,yMxMy)d'H n - 1 (x)dn n - 1 (y) 

Jm Jm 

+ 4 7 f {Vv E , VM )ip 2 dU n -\ 
Jm 

where ip = (X, vm) G H 1 (M) with J M ipdW 1 ^ 1 — by (|2.3p . Here Bqq stands for the second fundamental 
form of dSl, i.e., the Hessian of the distance function from the boundary and |-Ba/| 2 is the sum of the square 
of the principal curvatures of M. If <p is not continuous the boundary term J Mnaa Bqq,(vm, vm)*P 2 dV. n ~ 2 
is defined by the trace operator. We will return to the subject of second variation in section 4 and the 
role of the function (p. 

It is rather straightforward to show that if C 1 regular set E is a local minimizer of (|1.3j) . i.e., there 
is S > such that J(F) > J(E) for every \F&E\ < 5, then the second variation is positive semi-definite 

d 2 J{E)[p] > 

for every (p £ H 1 (dE) with J gE ipdW 1 ^ 1 = 0. The proof for this result is very similar to (pQ, Corollary 
3.4) and will be omitted here. As it was pointed out in the introduction, our main result Theorem 12.31 
deals with the question whether we may conclude the stability from the formula (|2.6[) . 

Our goal is to prove that a regular critical point, i.e. a critical point with C 1 -boundary, with positive 
second variation is indeed a isolated local minimum. 

Theorem 2.3. Suppose that E C with J^ue dx = m is a regular critical point of (|1.3[) which satisfies 

d 2 J(E)[ip] > 

for every ip € H 1 (dE) with J dE ipdH"^ 1 = 0. Then E is a strict local minimum and there are c > and 
6 > such that 

J(F) > J(E)+c„\FAE\ 2 
for every F € BV(£l) with JqUf dx = m and \FAE\ < S. 

In order to prove Theorem 12.31 we need to calculate the second variation at any regular set which 
satisfies the orthogonality condition, not just for critical ones. This is done in Proposition 14. II and thus 
we obtain the formula (|2.6|) as a corollary. 

It is well known that the functionals £ e T-converges in L 1 to £. As a corollary we obtain the following 
stability result. 

Corollary 2.4. Assume E is a regular critical point of J with positive second variation and denote 
U E = Xe — XQ\e- There is £0 and a family {u £ } £<£0 of isolated local minimizers of £ e with constraint 
J n u £ dx = m such that u E — > ue in L 1 (0) as e — > 0. 

3. Critical sets and Regularity of w-minimizers 

Let us briefly discuss about the regularity of local minimizers of (|1.3j) . Let us assume that we know, 
a priori, that our local minimizer is C 1 regular then we may use the Eulcr-Lagrangc equation (|2.4|) to 
deduce that it is in fact C 3,Q regular, as we pointed out in the previous section. Since our attention 
is in local minimizers, not in critical points, we will use the minimality itself to obtain the first order 
regularity. This will be done by using the theory of w-minimizers (of area). 
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Definition 3.1. A set E £ BV(Vl) is an w-minimizcr with constants A > and r > if for every G C f2 
with GAE C B r (x ) it holds 

P BAxa) (E) < P Br{xo) (G) + A\GAE\. 

We remark that the above definition differs slightly from the standard definition of w-minimizcr 
which appears in the literature e.g. in and [23) . where the measure of the symmetric difference 
|GA-E| is replaced by the volume of the ball \B r \. Definition 13.41 is clearly stronger and we may therefore 
apply known results from previously mentioned works. The motivation for the stronger definition is that 
for regular w-minimizcr we obtain a curvature bound. 

Proposition 3.2. Suppose that E G BV(£l) is an uj-minimizer with constants A > and r > and that 
E is C 1 regular. Then there exists a bounded function f such that E satisfies the equation 

(3.1) H M (x)=f(x) xeMnn 
in a weak sense. Moreover it holds < A. 

Proof. We may locally write M as a graph of a C 1 -function. Suppose that <p : D C W l ~ 1 — s- E is such a 
function and let rj € Cq (D) . From the w-minimizing property we obtain for every t it holds 

/ yJ\ + \V(f)\ 2 dz< [ y/l + \V^) + tVri\ 2 dz + At f \n\dz. 

JD JD JD 

Divide by t and let t — > to obtain 

f (V^Vt?) , . 

/ H dz < A 77 

v/i + iWF 

Arguing exactly as in Proposition 7.41 in [3] we conclude that E 1 satisfies (|3.1|) for some / with < 
A. □ 

We return to the C^-regularity of w-minimizcrs. The following regularity result in a smooth domain 
Q follows from a work by Griiter [llj . 

Theorem 3.3. Suppose that E is an uj-minimizer in f2. Then dE is C 1 regular outside a singular set T 
with Hausdorff dimension dim%(T) < n — 8 and M = dE D meets dfl orthogonally on dtt \ T. 

In [llj Griiter states the result for minimizers of the partitioning problem, i.e. p. 31) with 7 = 0. 
However he proves the result first by showing that the solutions of the partitioning problem are w- 
minimizers. Then he uses a result from |12j to prove the C^-regularity by using the fact that mean 
curvature is L p -integrable with p > n which follows from the corresponding Euler-Lagrange equation. 
Hence may apply this result to our w-minimizers which by Proposition ^. 2l have bounded mean curvature. 
The interior regularity follows from a work by Tamanini |23j . 

As in [Tj (Lemma 2.6) we have the Lipschitz continuity of the non-local part, i.e., for two measurable 
sets E, F C ft we have that 

(3.2) I [ \\7v F \ 2 dx- [ \S7v E \ 2 dx\ < C\FAE\ 

Jn Jn 

for some dimesional constant C. Here ve and vf are defined as in (|1.4|) . We may use an argument from 
[TTj to prove that local minimizers of (| 1 . 3|) arc w-minimizers in the sense of Definition 13.41 Although the 
proof is an exact copy from |llj we include it for the convenience of the reader. 

Proposition 3.4. Suppose that E G BV(Cl) is a local minimizer of (|1.3|) . Then it is an uj-minimizer. 
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Proof. Let G C il be such that GAE C B r (xo) for some r > 0. We may assume that \G\ < \E\ for in 
the case \G\ < \E\ wc argue similarly. Find a ball B r (yo), disjoint with B r (xo), such that 

\B r (y Q ) \E\ = \EH B r (x )\ -\Gn B r (x )\. 

By choosing G = GU B r (y ) we have that \G\ = \E\. The local minimality of E yields J(E) < J(G) and 
by (|3.2|) we deduce 

(3.3) Pba*«)( e ) + p ( E n Br(W))) < ^(.„)(G) + P(B r (lfo)) + C \GAE\. 

We may use the isoperimetric inequality to estimate 

P(EnB r (y )) > P(B r (y ))-C\B r (yo)\E\. 
Notice that \GAE\ < 2|GA£'| and \B r (y ) \ E\ < \GAE\. Hence, (JHSJ) implies 

P BAxo )(E) < P BAxo) (G) + A\GAE\. 

□ 

We may use the regularity theory for w-minimizcrs and the equation (|2.4I) to obtain that every local 
minimizer of ()1.3[) C 3,Q -regular outside a critical set T with dim^(r) < n — 8. This motivates us to define 
regular critical sets as critical sets with no singularities. 

It is well known that in higher dimension n > 8 the singular set even of a minimal surface might not 
be empty and therefore the regularity for w-minimizcrs is optimal. However by the flatness theorem of 
De Giorgi (see [10]) we obtain the full regularity if we know that the blow-up limit at every point on the 
minimal surface is flat enough. Since every blow-up limit of an w-minimizer converges to a minimal cone 
we obtain the following convergence result which concludes the section. The result can be obtained by 
following [2], [13] and [21] as explained in [24] with a few modifications. See also }23| . 

Theorem 3.5. Suppose that Ek are u> -minimizer s in f2 with uniform constants A and r. Suppose that 
E has C l ' a boundary in VL with a G (0, ^). // Ek — > E in L 1 then 

dE k -> dE in C 1 ' . 

In particular, dEk is C ,a -regular when k is large. 

4. Second variation formula 

In this section we calculate the second variational of the functional (|1.3[) . Our main task is to write 
the formula in a form where the quadratic structure appears. As in [I] we have to generalize the formulas 
from [7] and [22] and calculate the second variation at any regular set E, not necessarily critical, and 
with respect to any given flow. 

We assume that E C R" is C 2 regular by which we mean that M = dE n Vl is a C 2 -manifold with 
boundary. As in section 2 we consider a vector field X £ C 2 (f2;R") which satisfies the tangent condition 
(|2.2[) and the associated flow $:!lx (—t ,t ) —> Q, 

■ A 

-$ t =X($ t ), $(x,0)=x. 
ot 

We define the second variation at E with respect to the flow $ as 
where E t = <&{E,t). 
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A major technical challenge are the calculations for the local part of the energy Pq(-). We follow 
a slightly different path than ([4], Proposition 3.9) where the proof contains calculations for the general 
second variation formula of the perimeter in R n . Instead of differentiating the first variation formula 
in Proposition 12.11 we will first calculate the second derivative of Pn(Et) as it is done e.g. in [TU] and 
then use the divergence theorem in order to find the "right" formula. This saves us from differentiating 
the boundary term which appears in Proposition 12.11 Along the calculations it will become clear that 
we need to assume the orthogonality condition (|2.5[) on E in order to find the quadratic structure in its 
second variation. As we pointed out in the previous section the orthogonality condition appears in the 
regularity theory for aj-minimizcrs. It can thus be viewed as a part of regularity assumption on E. 

We recall the divergence theorem on M 

[ divA/AtfH"- 1 = f H{X,v M )dU n - 1 + [ (X,v*)dH n - 2 , 
Jm Jm JMndn 

where v* is the unit co- normal of M D <9f2, see Proposition 12. II Imposing the orthogonality condition on 

E implies v* = vq on M H dVl. We may thus write the integration by parts formula for / € C 1 (r2) as 

/ fdivMXdn n - 1 = -f(D T f,X)dn n - 1 f Hf(X,v M }dH n - 1 + [ f(X,v n )dH n - 2 , 
Jm Jm Jm Junen 

where D T is the tangential derivative on M , i.e., D T f = Df — (Df, Vm)vm- 
Proposition 4.1. Let E C K" be C 2 regular set which satisfies the orthogonality condition 



denote M = dED Q. Suppose that X satisfies the tangent condition (|2.2j) and denote X v = (X,vm)i / m 
and X T = X — X . The second variation of (|1.3[) at E with respect to X can be written as 

d2j(Et -Ln = [ (\D T (X,v M )\ 2 \B M \ 2 {X, VM ) 2 ) dU"- 1 - I . 

>M v JMndQ 

in— 1 / a injn— 1 / 



in— 2 



dt 2 l4 =° 

+ 87/ / G{x,y){X{x),v M ){X{y)^ M )dH n ~\x)dH n - l {y) 
Jm Jm 

+ 4 7 f (Vv E , v m ) (X, v M ) 2 dU n - Y - [ (H + A 1VE ) div T {X T ({X, v M ))) dU^ 1 
Jm Jm 

+ I (H + A 1 v E )&w(X){X,v M )d'U n - 1 . 
Jm 

Here Bqq stands for the second fundamental form of dQ, i.e., the Hessian of the distance function from 
the boundary and \Bm\ 2 is the sum of the square of the principal curvatures of M. 

Before the proof we remark that the above formula consists of the quadratic form (|2.6[) by defining 
ip := (X, vm) and two additional terms. It turns out that both of these extra will terms vanish when E 
is critical. Indeed, since the minimizing problem (|1.5p is volume-constrained it is natural to assume that 
the flow $ preserves the volume 

1^1 = \E\ 

where E t = <&{E,t). From the volume constraint we obtain the first order condition Jj|-Ei|| t _ = which 
yields 

div(X) dx= I (X, v M ) dU n - x = 0. 
Jm 



The second order condition Jp-|.Et||._ = yields 



/ div (div(A)A) dx = [ dW(X)(X, v M ) d 
Je Jm 
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These calculations are standard and can be found e.g. in [7]. Therefore if E is a regular critical set, 
i.e., H + Ajve is constant and the flow is volume preserving we obtain the quadratic form (|2.6|) with 
ip = (X,i>m)- Arguing exactly as in (pQ, Corollary 3.4) we deduce that every regular local minimizer of 
J satisfies d 2 J{E)[tp] > for every ip € C 2 {M) with ^tpdW 1 ' 1 = 0. 
We give now the proof of the Proposition 14. II 

Proof of the Proposition ^. 1\ Denote the acceleration vector field by Z = -§[z&t\ t _ = DX [X]. We treat 
the perimeter and the nonlocal part of the energy separately and denote 

E(t) = Paa(Et), F(t) = 7 / \Vv Et \ 2 dx. 

Jn 

For simplicity we write v = vm- 

The value for F"(0) is calculated in [7j in the periodic setting. However, as it was pointed out in 
Remark 2.8 in [7], the Neumann boundary case does not produce any new terms to it and the formula 
for F"(0) is the same as in the periodic case. Recalling the condition (|2.2|) we have from [7] 

F"(0) = 8 7 / / G{x,y){X{x),v){X{y) 1 u)dU n - 1 {x)dH n - 1 {y) 
^\ Jm J m 

+ 4 7 / 61v{v E X){X,u)dH n ~ 1 . 



JM 

The value of E"(0) is well known in the following form, [10] (Theorem 10.4), 

E"(0) = f (div A/ Z + (div M X) 2 + \{D T X) T v M \ 2 - Ti-(D T X) 2 ) dH n ~\ 
Jm 

where divm X = div X — (DXv,v) and (D T X)ij = (D T Xj,ei) = DX — (DXv) <E> v. By the divergence 
theorem we have 

(4.2) E"(0)= f \(D T X) T v\ 2 + (div M X) 2 -TY((D T X) 2 )+H(Z,v)dH n - 1 + f (Z, m )dH n ~ 2 , 
Jm JMnn 

where we also used the orthogonality to conclude v* = Uq on M n 90. Notice that this is the second 

variation formula for the perimeter Pn(E t ) given by any vector field X. 

We have to manipulate (|4.2[) in order to find the quadratic structure in it. We proceed by writing 

X = X v + X T where 

X v = {X,v)v. 

The goal is to decompose every term in (|4.2[) by using X = X v + X T . We will first treat the term 
J M \(D T X) T is\ 2 dH"' 1 in Since Dfi/ = 0wc have {D T X v ) T v = D T {X,v) and therefore 

\{D T X) T V \ 2 = \(D T X v ) T v\ 2 + 2{{D T X U ) T V ) ■ {{D T X T ) T v) + \{D T X T ) T v\ 2 
^ 4 ' 3) = \D T {X,v)\ 2 +2D T {X,v) ■ ((D T X T ) T v) + \{D T X T fv\ 2 . 

Since divM X v = H(X, v) the second term J M (divM X) 2 dW 1 ^ 1 can be written as 
(div M X) 2 = (div M X v + div M X T ) 2 

(4.4) 

= H 2 (X, v) 2 + 2H(X, v) div M X T + (div M X T ) 2 . 

For the third term we will use the equalities (D T X U ) 2 = (X, v) 2 (Dv) 2 + (D T (X, v) ■ v) v ® D T {X, v) and 
(D T X T )(D T X V ) = (X, v)Dv D T X T . Hence, we deduce 

Tr(D T X) 2 = Tr(D T X u ) 2 + 2Tr {{D T X T ){D T X U )) + Ti(D T X T ) 2 

^ = {X,v) 2 \B M \ 2 +2(X,u)Tr(DuD r X T ) + Ti(D T X T ) 2 . 



in 
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We treat the fourth term J M H(Z, v) dH n 1 by writing (here we denote Z T = D(X T )(X T )) 

(Z, v) = (DXX, v) = ((DX(X U ) + D(X„)(X T ) + D(X T )(X T )) , v) 
= (X, v)(DX v, v) + D T (X, v)-X T + (Z T , v) 

(4.6) 

= (X, v) divX - (X, v) div M X + D T (X, v)-X T + (Z T , v) 

= (X,v) divX- dw M {X T (X, v)) - H(X, vf + 2X T ■ D T (X,v) + (Z T ,v), 

where in the last equality we have used divjv/ X v = H(X, v). 

For the last term J M Q (Z,i>n) dH n ~ 2 we notice that the tangent condition (|2.2[) on X implies 
(X, uq) = on M n dSl. In particular, D(X, vq) ■I=0onMn dfl. Hence the last term becomes 

(4.7) f (Z, v Q ) dH n - 2 = f (DXX, u Q ) dH n -' 2 = - f (Dv n X, X) dH"- 2 . 
JMndn JMndn JMndn 

From now on we will use the notation (DvqX, X) = Bqq(X, X). 
We use flU), (|P|) . (gU) and (gTTJ) to rewrite flO) as 

E"(0)= [ \D T (X,v)\ 2 -\B M \ 2 (X,v) 2 dU n - 1 - [ B d n(X,X)dH n - 2 

JM JMndn 

+ 2 f D T (X,v) ■ ((D T X T ) T v)dU n - 1 +2 f H ((X, is) div M X r ) dH 71 - 1 

JM JM 

(4.8) -2 / {X,v)Tr (Dv D T X T ) dU n ~ x +2 / H (X T ■ D T (X, u)) dH^ 1 

JM JM 

+ I H {2X T ■ D T (X, v) + (X, v) divX - div M {X T (X, u))) dU^ 1 

JM 

+ [ (div M X T f + \(D T X T fv\ 2 -Tr{D T X T ) 2 + H(Z T ,v) dU n ~ x , 

JM 

where Z T = D(X T )(X T ). 

We continue by treating the mixed terms in (|4.8[) . We integrate by parts the third term in (|4.8[) and 
deduce 



2 I D T (X,v) ■ ((D T X T ) T v)dW 

JM 



(4.9) =-2 / {X,v)dw M ({D T X T ) T V )dU n - 1 +2 / (X,v)((D T X T fv,v Q ) dH n ~ 2 

Jm JMndn 

= -2 [ {X,v)div M ({D T X T ) T v)d'H n - 1 -2 [ {X, v) (Dv X T , v n ) dH n ~ 2 , 
jm JMndn 

where we have used = D T {X T , v) ■ = (D T X T i/q, v) + (Dv X T , vq). 

For the fourth term in (|4.8p we observe that, by the orthogonality condition, X T vanishes on Mndfl 
and therefore integration by parts yields 

(4.10) 

2 f H(X, v) div M X T dU n - x = -2 [ D T {H(X, v)) ■ X T dH"- 1 +2 [ H(X, u)(X T ,m) dU n - 2 
jm Jm JMndn 

= -2 I HX T - D T (X,v)dH n ~ 1 -2 (X, v)D T H ■ X T dU n ~ x . 

JM JM 

Next we introduce a notation 5i for partial derivative on M, i.e., Sig = d Xi g — (Dg,v)vi for any 
smooth function g. We use the facts SjUi = SiVj and Y^h j{^3^i v i)-^-T = ^Zi j(^i p j)-^T ( Lemma 10.7 in 
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[TO] ) to obtain 

div M {{D T X T ) T v) + D T H ■ X T + Tr (Dv D T X T ) 

n 

= »i) + (s^xi + 6^x1 

(4.11) 



id 
n 

= J2 hHXi V j) = A M (X T ,v) = 0, 

i,j 

where Am = X)"=i IS ^ ne Laplacian on M . 

Together with P~9|) . (|4TT))) and (|47TT|) the formula (|4~8)) becomes 



E"(0)= f \D r (X,v)\ 2 -\B M \ 2 (X,v) 2 dU n ~ l 

JM 

- [ B dn (X,X)dH n - 2 -2 f (X,v)(DvX T ,v Q )dU r 
JMndn JMndn 



(4.12) 

+ / H((X,v)divX -div M (X T (X,v))) dU" 1 - 1 



M 



[ (div M X T ) 2 + \(D T X T ) T v\ 2 -Tr{D T X T ) 2 + H (Z T , v) dW 

JM 



I M 

Let us treat the last row in the previous formula. We claim that 



(4.13) / (div M X T ) 2 + \(D r X r ) T v\ 2 -TY{D T X T ) 2 + H(Z T ,is)dH n - 1 + f (Z T ,v n ) dH n - 2 = 0. 
JM JMndn 

Indeed, the formula (|4. 13[) is nothing but the second variation formula (formula (|4.2[l ) of the perimeter 
with respect to the vector field X T . As we noticed before X T vanishes on M n dfl and therefore the flow 
associated by X T leaves M unchanged. Hence, we have (14. 13[) . 

Since X T vanishes on M n Oft, we have D(X T , vq) ■ X T = on. This yields 

(4.14) / (Z T ,v n )dH n - 2 = [ (D(X T )X Tim )dH n - 2 = - f B dn (X T , X T ) dW 1 ' 2 . 
JMndn JMndn JMndn 

Moreover, since (v,Vq) = on M n dtt we get = D(v,vq) ■ X T = (Dv X T) vq) + (DvqX t ,v) and 

therefore 

(4.15) 

2/ (X,v)(DvX T ,v u )dU n - 2 = -2 [ (Dv n X T ,X v )dU n - 2 = -2 [ B dn (X T ,X„)dH n ~ 2 . 
JMndn JMndn JMndn 

We combine (|4.1j) and (|4. 12[) . together with (I4.13[) . (|4.14l) and (|4.15|) and write the second variation 

as 

-/ (\D T (X,v u )\*-\B M \ 2 (X,u M ) 2 ) dH*- 1 - f B &n {v,p)(X,p) 2 dH n ~ 2 
J m ' JMndn 

I H({X,v)divX-dW M (X T {X,v))) dH 71 - 1 

JM 



dt 



2 li=0 



(4.16) 

+ 87/ / G(x,y)(X(x),v)(X(y),v)dH n - 1 (x)dH n - 1 (y) 

JM JM 

+ 4 7 f dW(v E X) (X, v) dU n - 1 . 

JM 
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As in [T] we integrate by parts to obtain 

dW(v E X)(X,v)dH n - 1 = [ (Dv E ,X)(X,v)dH n - 1 + [ v E div(X) (X, v) dW 

M JM JM 



[ (D T v E ,X T )(X > u)dH n - 1 + [ (Vv E ,v)(X,v) 2 dH n - 1 + [ v E divpQ (X, v) dH r 

JM JM JM 



in—1 

v E uiv^yi ) v j ur 

' M JM JM 

in— I i / /Y7„._ ,,\/v ,,\2 Jun— 1 i / A-„,t v\ I v ,.\ jun—1 



v E div(X T (X,v))dU n - 1 + / (Vv E ,iy)(X,iy) 2 dH n - 1 + / v E div(X) (X, v) dU n -\ 

M JM JM 

where the last equality follows from the fact that X T vanishes on MOdil. Combining this to (|4. 16[) yields 
the result. □ 

5. VF 2 ' p -MINIMALrTY 

In this section we prove that a regular critical set E with a positive second variation is a strict local 
minimizer among sets which are regular and close to E in a strong sense (in W 2,p topology). This result 
is stated in Proposition 15.21 and is interesting in itself. The idea is to take an arbitrary set F near E and 
to construct a volume preserving flow $ such that &(E, 1) = F. We use Proposition 14. 1 1 to estimate the 
second variation at E t = $>(E,t) at any time t <G [0, 1]. Notice that by the assumption on E we know 
that the second derivative of t i-> J(E t ) at t = is strictly positive. We then use the fact that E t are 
close to E in W 2,p topology to deduce that the function t H> J(E t ) is in fact strictly convex. We note 
that in order to use Proposition 14.11 all the sets have to satisfy the orthogonality condition. 

We begin by a simple compactness argument. The proof is exactly the same as Lemma 3.6 in pQ 
and will therefore be omitted. 

Lemma 5.1. Suppose that E is a critical point with positive second variation. There is cq > such that 

d 2 J(E)[<p] > co\\<p\\ 2 H i( dE) 
for every ip g H 1 (dE) with § dE LpdU n - 1 = 0. 

We define W 2,p and C 2 distances between regular sets E, F C as 

||£',F||w2,p := inf{||$ - Id\\ W 2, P ^ \ $ : E ->• F, C 2 - diffeomorphism } 
\\E,F\\ C 2 := inf{||$ - Id\\c*(fi) I $ : E -> F, C 2 - diffeomorphism }. 
Here is the main result of the section. 

Proposition 5.2. Let E be as in Theorem \2.3\ and p > n. There is 8 > and a constant c\ > such 
that for every Fc!! with \F\ = \E\ which satisfies the orthogonality condition and \\E, F\\ W 2, P < 5 it 
holds 

J(F) > J(E) + Cl \FAE\ 2 . 



We first prove some technical lemmata and then give the proof of Proposition 15.21 at the end of the 
section. Most important of the lemmata is Lemma 15.31 where we construct a vector field X and a flow $ 
which takes the given set E to a regular set F nearby. The proof is technical and long and it is therefore 
divided into three steps. The difficulty lies in the fact that the flow needs to satisfy both orthogonality 
and the volume constraint. We divide the proof by first constructing a primitive flow which satisfies the 
orthogonality condition and then, at the third step, we correct it so that finally also the volume constraint 
is satisfied. 



MINIMALITY VIA SECOND VARIATION FOR MICROPHASE 



13 



Lemma 5.3. Let E be as in Theorem \2.S\ Suppose F is C 2 regular set which satisfies the volume 
constraint \F\ = \E\, the orthogonality condition (|2.5|) and \\E, F\\yy2, P < 5, where p > n. When S is 
small enough there is C 2 regular vector field X with uniformly bounded W 2 ^ '-norm and associated flow 
$ 

^*(x,t) = X{$(x,t)), $(x,0)=x, 

such that at every t € [0,1] the set E t = <f>(E,t) satisfies the volume constraint \E t \ = \E\, and the 
orthogonality condition (|2.5p . Moreover the flow takes E to F, i.e., E\ = F . 

Proof. Step 1: First we construct a primitive vector field Z which has the property that the associated 
flow $ z : O x (-t , t ) -> O, 

d 

— §z = Z($ z ) with $z{x,0) = x 

is such that, at any t e (—to, to), the set &z(E,t) satisfies the orthogonality condition. We will build the 
vector field only in a neighbourhood of the relative boundary M = dE n fi, which we denote by U. By 
neighbourhood of M we mean a connected set U which contains M and is open with respect to Ct. We 
begin by defining a function which plays the role of a signed distance function of the set E with respect 
to the domain Q. 

Since E is C 3 ' Q regular and satisfies the orthogonality condition there is a C 3,a regular function 
dn.E '■ U M which satisfies 

(i) dn,E = on M and d^,E < in E, 

(ii) \Vd n , E \ >O0inf7, 

(iii) ^Oonfflnt/. 

Such a function can be found e.g. by minimizing the energy 

/ \Du\ 2 dx ueW 1 ' 2 (U\E) 

JU\E 

with respect to Dirichlet boundary condition u = on M and u = 1 on (dU nti)\E. Here we assume 
without loss of generality that the boundary dU n il meets dft orthogonally. We may choose dn.E m 
U \ E to be the minimizer of the above problem. One may then extend this function naturally to whole 
U. 

It follows from the above conditions (i) — (iii) that for small values of t € R the level sets 

{xeQ \ d QtE (x) = t} 

form a foliation of U with C 3,a regular sets which meet dd orthogonally by the Neumann boundary 
condition (iii). We will define the primitive flow $^ such that 

$z(E,t) = {x g A | d u>E (x) < t}. 

For future purposes we need yet to construct the trajectories of the flow. 

To that aim we construct a vector field Z' which defines the trajectories. We begin by defining 

Zi := Vdn,E in U. 

To make sure that we will not loose the orthogonality in the forthcoming steps we correct Z\ near the 
boundary dfl. This will be done by "copying" the geometry of dil into its neighbourhood. We denote 
the neighbourhood of d£l by U'. 
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Since <9f2 is regular we may define the standard projection to <9f2, II : U' — > dCl, as II(x) = y x , where 
y x £ dQ is the unique point for which 

d(x,fi) = inf \x-y\ = \x- y x \. 
yen 

Every point x £ U' can therefore be written as x = n(x) — d(x,Q)vQ(H(x)) where vn{y) is the outer 
normal of il at y £ dil. We define the field Z2 such that Z2 — Z\ on the boundary 9fi l~l U and 

Z a (iO = Zi(II(|/)) for y e t/nt/'. 

The advantage of Z2 is the following geometrical fact. If we know that the points x, y £ 9f2 n U lie on 
the same trajectory, then for small h > also the points x — hvpi(x) and y — /ifn(y) he on the same 
trajectory. We will use this fact frequently. 
Finally we define Z' as 

Z' = r]Z x + (f - T))Z 2 

where 77 £ Co°(f2) is a standard cut-off function such that 77 = near 9fi and 77 = 1 outside [/'. By 
choosing the cut-off function properly we can make sure that Z' is close to Vdo^, i.e., 

\\Z' - Vdn,B|U»( A /) < e. 

For every point x £ M we define the primitive flow x(t) such that the trajectories are given by Z' 
and x(t) is the unique point on the intersection of the level set {x £ fl \ (1q,e(x) = t} and the trajectory 
passing through x. In other words x(-) is such that x(0) = x and 

(5.2) dn, E (x(t))=t. 

We denote the flow by $z(x, rj) and the associated vector field by Z. Since Z' is close to Vdn,_B on 
M = {x I dn,_e(x) = 0} we have that 

(5.3) \{Z{x),v M )\>{\-e)\Z{x)\ for all x £ M. 

The flow $z satisfies the orthogonality condition by construction. In fact the construction implies 
that every point y £ d£l n U can be written as y = $z(x, t) for some t £ (—t 0} to) and it holds 

(5.4) D<S> z (x,t)v Q (x)=v Q (y). 

This is a priori a stronger than the orthogonality condition since it implies that the vector D<&z{x, t)i/a(x) 
is the unit co-outer normal of the set <&z(M,t) and it equals Vq($>z(%, t))- 

To verify (|5.4j) write y = $^(x,t) and notice that by the construction, for every small h > 0, the 
points x — hvn(x) and y — hvn(y) lie on the same trajectory. Therefore there is t% such that 

y - hv n {y) = $ z (x - hv n (x),t h ). 

In particular, we deduce from (|5.2j) that da t E(y — hvo,{y)) = th- 

Notice that dn,B(y) = da,E(^z(x, t)) =7;. Therefore by the orthogonality of the set {x | c^o,b(x) = t} 
it holds 

rfn,£(y - hu a (y)) = t + o(h). 
Consequently th —t + o(h) and (|5.4j) follows. 



Step 2: Next we change Z to Y so that the new flow, denoted by $>y, satisfies $y(-E, 1) = F and 
the orthogonality is preserved. 
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In the neighbourhood U of M the map &z\Mx(-t ,t ) — ► U is C 3:Q -diffcomorphism and every point 
y £ U can be written as y = &z(x,t) for some x £ M and for t = dn t E(y)- We may therefore define 
implicitly a projection tt : U M such that 

®z(ir(y),<ki,E(y)) = v- 

Since Z and dn,E are C 3,Q regular, also 7r is C 3,a regular. 

By the assumption on F there is C 2 -diffeomorphism : E 1 — > F with ||\& — 7d||iy2.p < S. Without 
further mention we always assume S > to be small. We define a map S : M — > M as 

= vr(*(x)). 

The tangential differential of S is 

A-S^x) = A-7r(*(a)) £> r *(x) on x £ M. 

From the regularity of 7r and from the fact 7r(x) = x for x <G M we conclude that |D T 7r('J'(x))| > c. Since 
||D T v E'(x) — /||i<» < S we conclude that S is C 2 -diffeomorphism and the bound ||4 , || H a2, p < C implies 

(5.5) WS^Wwhm) <C. 

We define further a projection on F Mp := dF P\ £1, np : U ^ Mp, as 

n F (y) = y(S- 1 (TT(y))) 

which labels for every point y £ U a unique point z y on Mj? on the same trajectory. Trivially the map 
TtF is constant along the trajectories of <£>z and the value Tp{y) := dn.E(^F(y)) denotes the time needed 
from M to Mp along the trajectory passing through y. We define the vector field Y by 

(5-6) Y(y)=T F (y)Z(y) 

and denote the associated flow by $y. We note that we may write 3>y explicitly as 

(5.7) $y(x,t) = $ z (x,T F (x)i) for x € M. 

We denote by F t the set enclosed by $y(M, t). Since Tp(x) is the time needed from x £ M to Mp 
it follows that Mp = {$y(x, 1) | x G M} and therefore F\ — F. Moreover we have that ||Tf||w3,i> < 
\\F,E\\ W 2, V <5. 

We have to make sure that at any time t £ (0, 1) the set F t satisfies the orthogonality condition. As 
in the previous step we will show that for any y £ $y (M, t) n 90 it holds 

(5.8) D9 Y (x,t)va(x) = va(y) 

where y = <3>y(x,i) for some x £ M. This implies the orthogonality. 
To show (|5.8j) we use (|5.7j) to calculate 

9$ 7 

£>$y (x, t)va (x) =D*z(x, T F (x)t) v n (x) + i-^- (x, T F (x) (x)t) ( VT F (x) , m (x)) . 
Hence by (|5.4[) we need to show 

(5.9) (VT F (x),i/ n (x)) =0 for x e ill n 90. 
Since Tp(y) = dn^E(Trp(y)), (|5.9j) can be written as 

(Vd aB (7r F (x)),i}7r F (xH(2;)) = for x £ M n 90. 
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We claim that Dirp{x)vo,{x) = vq_(tip(x)). Then (|5.9j) follows from the fact that da,E satisfies the 
Neumann boundary condition 

ddn,E _ Q 

To show D-kf{x)vq,{x) = vn(y) where y = 7Tp(x), recall that irp is a map which labels for every x' € U 
a unique point y' € Mp on the same trajectory. Therefore the points x — hva(x), y — hv^{y) and 
ttf(x — hi>n(x)) all lie on the same trajectory. It follows from the orthogonality of F that 

tt f (x - hv n (x)) = y- hv n (y) + o(h) 

and the claim follows. 

At the end of the step we will write down some useful properties of the flow $y and the field Y for 
the next step. First of all, it follows from the definition of Y (|5.6[) . (|5.3I) . (|5.5|) and the fact that Z is C 2 
regular that 

(5.10) OTfI < \Y\ < C Tp and \\Y\\ W 2, P < C\\T F \\ W 2. P < 5. 
In particular since Tp is constant along the trajectories (|5.10[) yields 

(5.11) C _1 |T F (s)| < \Y(&r(x,t))\ < C\T F {x)\ for all x e M 
and we may also estimate the divergence by 

(5.12) |div(y)| = |!Z>div(Z)| < C\Y\. 
Since $y solves 

d 

— $y (x, t) = Y($ Y (x, t)), and $ r (x, 0) = a; 

we may estimate, by differentiating the above equation for $y twice and by using the second inequality 
in (pnO]) . that 

||$y(-,t) - i"d|| W 2, P(M) < C||T F || W 2. P < 5. 
The above estimate together with (|5.3[) imply 

(5.13) |(y,^ Ft )|>(l-2 £ )|F|. 

Step 3: We modify Y to X so that the new flow will finally be volume preserving. We note that 
the sets 

F t = $y (E,t) 

satisfy everything except the volume constraint but the final set = |F| has the right volume. The 
idea is to study the evolution of F t \ E and E\F t . If the volumes of these two sets do not match we give 
the other one a little boost. More precisely we would like to find for every t <G [0, 1] a value s(t) such that 

\F t+s{t) \E\ = \E\F t \. 

The problem is that this would break the regularity of the flow. 

To overcome this problem we choose a smooth function / defined on M such that that / ~ XMnF- 
To be more precise, we choose / such that it satisfies 

(i) / = 1 on M n F and / < 1 on M \ F, 

(ii) / = on M \ F e where F £ := {x 6 il | dist(x, F) < e}, 

(iii) (V/, v a ) = on M n 90. 
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For every t £ [0, 1] we define a flow $> t : M x [— t, 1 — i] — > f2 such that 

$t(x,s) = $y(x,t + s/(x)) 

and the set enclosed by <E» t (M, s) by F tjS . We denote the associated vector field by X t which we may 
write as 

X t (y) = f(ir(y))Y(y) 

where 7r is the projection on M defined in the beginning of step 2. We note that Ft,o = F t is the set 
given by the flow $y . Moreover since / = 1 on M H F we have 

(5.14) $ t (x,s)=& Y (x,t + s) for xeA/nf. 

The idea is to correct the volume of F t such that the new set, denoted by E t , will satisfy the volume 
constraint. To that aim we define a function 

p(t,«) :=\F t ,.\-\E\. 

The correction value s(t) is then implicitly defined as tp(t, s(t)) = and E t := F t s ( t y To be more rigorous 
for every t £ [0, 1] we define s(t) as 

(5.15) s(t) := sup{s G [-t, 1 - t] : <p(t, s) < 0}. 

First we have to make sure that s(t) is well defined. 

Firstly we show that for every t £ (0, 1) the function s >->• ip(t, s) is strictly increasing, i.e., 

(5.16) ^ > 0. 

as 

From (|5.1ip and (|5.13[) we conclude that the growth of the volume \Ft \ E\ do not degenerate 

(5.17) ±\ Ft \E\= [ {Y,v Ft )dU n ~ x ~ / \T F {x)\dU n -\ 

at JdF t \E JMCiF 

By notation a ~ /3 we mean that C _1 a < /3 < Ca for some constant C. In particular the above estimate 
implies 

(5.18) \FAE\ =2\F \E\ =2 f^-\F t \E\dt~ f \T F (x)\ dU n ~ x . 

Jo at J MnF 

Similarly 

(5.19) ^-\E\F t \= f \(Y, vf,) \ dT-C 1 ^ 1 ~ f \T F (x)\ dH^ 1 - \FAE\. 

at JdF t r\E J M\F 

Since the vector field X t has the same regularity as Y we may estimate the growth of \E\ F t s \ for 
every t £ (0, 1) and every s £ (— t, 1 — t) exactly as we with (|5.17[) and obtain 



|-|£7\f ti ,| = - / (X,u Ftt .)dH n - 1 <C [ f(x)\Y(x)\dh n - 1 <C [ f{x)\T F {x)\dW 
US JdFt „nE Jm\f Jm\f 



ldF tiS nE Jm\f Jm\f 

We note that T F (x) = for x £ M D <9-F. It is therefore clear that the above inequality and the condition 
(ii) on / imply 

(5.20) ^- s \E\F t , s \<Ce\FAE\. 

To prove J5T6]) we write ip(t, s) = |Ft lS \J5|-|^\Fi, a |. From ([5T4]) it follows that \F t , s \E\ = |F t+a \J5|. 
From (|5T7) . (|5~18f and ([5^20)) we conclude that 

*) = £ (|F i+s \ 25| - |25\ F t , a |) > (c - <7 £ )|FA£|. 
as as 
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Therefore (|5.16|) follows by choosing e to be small enough. 

Notice that s(0) = s(l) = 0. This follows simply from the fact \F\ = \E\ and (|5.16|l . To be sure that 
for every t the value s(t) indeed exists we claim that for t G (0,1) we have 

(5.21) -t<s(t)<l-t. 

To that aim we again write ip(t, s) = \F t+s \ E\ — \E \ F t . s \- Since Fq = E and Fi = F we conclude that 

\F tt _ t \E\=0 and \F t>1 _ t \ E\ = \F \ E\. 
Similarly since f(x) < 1 on M \ F it holds that 

\F t ^ t n E\ > and |F M _ t n E\ < \F \ E\. 

Hence ip(t, —t) < and tp(t 7 1 — t) > 0, and (|5.21[) follows. Since the flow $y is uniformly W 2 ' p regular, 
the function ip is C 2 regular. Therefore the condition (|5.16|) and implicit function theory imply that 
(|5.15p defines a function s(-) which is C 2 regular. 

In order to conclude that t + s{i) defines a paramctrization of the interval [0, 1] we show that for 
every t g [0, 1] 

(5.22) s' (t) > -1 + fj,, 

for some /i > 0. Similarly as we obtained f|5. 19[) we use the regularity of X t to deduce that for every 
t G (0, 1) and s G (— t, 1 — t) we have 

l\E\F t , a \= [ KY^FtJldH"- 1 >c f ^{x^dW 1 - 1 >c\FAE\. 

M JdF t , s r\E J M\F 

The above estimate, (|5.17|) and (|5 . 20[) yield 

= ±<p(t, s(t)) = j f (\F t+s{t) \E\-\E\ F tAt) |) 

< (Cmax{l + s'(t), 0} - c + Csup \s'{t)\e) \FAE\. 

Consequently we obtain (|5.22p when e > is small enough. 
We finally define the flow $ : X [0, 1] -> as 

(5.23) t) = $ t (x, s(t)) = 3> Y {x, t + s(t)f(x)) for x G M, 
and extend it smoothly to fi. The sets E t are defined as 

(5.24) E t :=$(E,t) = F tta(jt) . 

By construction Et satisfies the volume constraint at any time t G [0, 1]. It also satisfies the orthogonality 
condition, which can be checked similarly as it was done in (|5.9[) by using the Neumann boundary 
condition (V/, i^n) = on M n dQ. 

The map 3>mx[o,i] ~> FAE is C 2 diffeomorphism and we may define projections it : FAE —> M and 
T : FAE ->[0,i\ such that for every y G FAE 

<P(n(y),T(y))=y. 

It is clear that the projections are C 2 regular with uniformly bounded W^' p -norms. The vector field X 
can be written on FAE as 

(5-25) X(y) = (1 + s'(T(y)) f(n(y))) Y(y). 
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It is then clear that X is C 2 regular with uniformly bounded W 2 ' p -norm. □ 

In the next lemma we study the regularity properties of the flow $ and the vector field X constructed 
in the previous lemma. The reader is encouraged to skip it in the first reading. 



Lemma 5.4. Suppose that E, F, $ and X are as in Lemraa \5.S\ and denote M t := dE t H fi. The vector 
field X € C 2 (fl;R n ) satisfies 

(5.26) \FAE\<C [ \{X,v M )\clM n - 1 

Jm 

and 

(5.27) \\X Tt \\ m(Mt) <C\ \(X,u t )\ \ m{Mt) , \\dW(X)\\ L2{Mt) < C\\(X,v t )\\ LHMt) 
for some constant C . 

Proof. We begin by making some remarks on the regularity of the flow $ and the vector field X. Notice 
first that by (|5.25|) the vector field X can be written as 

X(y) = <j>(y)Y(y) 

where <j> is a C 2 regular function with uniformly bounded W 2,p norm and C^ 1 < <j> < C ' . The estimate 
\Y(y)\ < C\T F (y)\ proved in ([5TT0|) yields 

C-^Tp (y)\ < \X(y)\ < C\T F (y)\ for all y e FAE. 

We recall that for x € M, Tf{x) denotes the time needed from x to Mf along the primitive vector field 
Z and HTfH^.p < C\\F, E\\ W 2, P . Since Tp is constant along the trajectories it follows from the above 
estimate that for every t e [0, 1] 

(5.28) C _1 |X(a;)| < \X(®(x,t))\ < C\X(x)\ for all x G M. 

Denote <&*(•) = $(•, t) and by vu t the normal vector of Mt. For every t S [0, 1] the map $(•, t) is C 2 
diffcomorphism. We may show, exactly as we did (|5 . 1 3[) in Lemma |5.3[ that 

(5.29) \\*(;t)-Id\\ W 2, HAI) <C\\X\\ W 2, P <C\\T F \\ W 2. P 

and 

(5.30) \{X,v Mt )\>(l-3e)\X\. 

We proceed by showing (|5.26p . We use the same argument as we did with (|5.18[) and (|5 . 19[) and 
calculate 

d , „ . „, d 
~dV 

where we have used the regularity of the flow, i.e., | Ja/$*| < O, (|5.28p and (|5.30p . The previous calculation 
then yields 

I FAE I = 



\E t AE\ = ^- (\E t \E\ + \E\E t \)= f \{X,u Mt )\ dW"" 1 < C [ \(X,v M )\dH n ~- 
dt jM t Jm 



f ^-\E t AE\dt<C f \{X,u M )\dn n -\ 
Jo dt Jm 

The first estimate in ()5.27[) follows from (|5.29p and (|5.30[) and from the VF 2 ' p -regularity of the sets 
Et- In fact the proof is exatcly the same as the one of Lemma 7.1 in [I] and it will therefore be omitted 
here. 

To prove the second inequality in (|5.27l) we estimate the divergence by using (|5. 12[) 

Idivpoi = \dw(<j>Y)\ < |(v0,r)| + |0div(y)| < |V0||r| + c\cj>Y\ < c\x\. 
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Hence, the second inequality in (|5.27|) follows from (|5.30j) . □ 

In the following lemma we show by continuity that in a neighbourhood of a critical point with 
positive second variation the quadratic form remains positive (|2.6j) . Notice that we have to compute the 
quadratic form for sets which are not critical. 

Lemma 5.5. Suppose that E is as in Theorem 1 2. 3\ and p > n. There is S > such that for any C 2 -regular 
set F with \\E, F\\ W 2, P < S and every ip G H 1 (Mp) with J M (pd'H n ~ 1 = it holds 

d 2 J(FM > ||M& 1(aF) . 



Here Mf = dF n f2, d 2 J(F)[ip] is defined in (|2.6p and cq is the constant from Lemma \5.1[ 

Proof. We argue by contradiction and assume that there are F k with \\E, F\\w2, P — > and ip k g H 1 (M k ) 
with J Mk tp k d'H n ~ 1 = 0, which by scaling we may assume to satisfy ||<Pfc||if 1 (A/ fc ) — 1: such that 

(5.31) lim d 2 J(F k )[p k ] < ^ 

A;— >oo z 



where M k = dFk D £1. We recall the formula for d J(F^ 



d 2 J{F k )[ Vk ] = f (\D Tk <p k \ 2 - \B Mh \ 2 <p 2 ) clU"- 1 - f B an {v k ,v k )<pldU n - : 
JM k JM k ndn 

n—l /_-\ Ji_/n— 1 / 



(5.32) +87/ / G(x,y)<p k (x)p k (y)dH n - L (x)dH n - L (y) 



2 jn/n-1 



+ 47/ (Vv Fk ,v k )vlm 

where v k is the unit normal of F k . We will show that there is cp € H 1 (M) with H^Hjyirjvf) = 1 and 
J^jipdW 1 - 1 = such that 

lim d 2 J{F k )[p k \ > d 2 J(E)[ V ]. 

k— >oo 

up to a subsequence. This and (|5.31[) then contradicts Lemma l5.il 

Since \\E, F\\-^2, P — > there arc C 2 -diffcomorphisms $ k : E — > F k such that — Zd| | w 2 .^ — ^ 0. 
By compactness there exists ip e H 1 (dE) with Hv'Hhuob) = 1 and J gE (p dW 1 ^ 1 = such that, up to a 
subsequence, 

fk ° &k f weakly in H {M), 



where M = dE fl f2. In particular <^fc o <I> fc — > tp strongly in L (M). We also conclude that v k o $ fc — > v 
uniformly on M where v is the unit normal of E. Therefore by the weak lower scmicontinuity we obtain 
the convergence of the first term in (15.321) 



lim / \D Tk p k \ 2 dU n - x > [ \D T tp\ 2 dH r ' 

k ^°° JM k JM 



_/n—l 

We proceed by observing the following convergences 

B Mk o$ k ^B M in L p (M) and v Fk -> v E in C ha (Q). 



Indeed, the first one follows trivially from the iy 2,p -convergence of F k and the second one follows from 
the uniform C 1,Q -regularity given by the equation (|1.4[) . Therefore we obtain the convergence of the 
second and the last term in (|5.32[) 



/ \B Mk \ 2 p\ dU n - 1 -> / \B M \ 2 p 2 dU n -- 

JM k JM 
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and 



J M k Jm 
By compactness of the trace operator we have that <p k ° &k — ► <P on L 2 (M n (957). Hence we deduce 
the convergence of the boundary terms 



JM k ndn JMndQ. 



We have yet to show 

(5.33) / / G(x, y)tpk(x)ifk{y) dH n ~ 1 (x)dH n ~ 1 (y) -> / / G^y) V {x)^y)d^- x {x)dH n -\y) 

JM k JM k Jm jm 

to conclude the proof. Denote by Jm^r the Jacobian of $fc on M and notice that Jm^u — > 1 uniformly 
on M. For simplicity we denote ifk = <fi k ^fc- Firstly we deduce 

G(x,y) i p k (x)cp k (y)dH n - 1 (x)dH n - 1 (y) - / / G(x,$ fc ( 2 /))^(x)^ fc (y)JM$ fc (y)dH"- 1 (x)dH"- 1 (y)| 

M k JM k JMJM 

< I \JM*k(v)<Pk{y)\ I \G^ k {x)^ k {y))^ k {x)J M ^ k {x)-G{x,^ k {y)) V {x)\dH n ~\x)dH n -\y) 

JM JM 

/ 2 \ 1 / 2 



< C 



\ 1/2 / , x 1/2 

$fc(i/))Jjtf*fc(a:) - G(x, $ fc (y))| 2 dU^dW 1 - 1 ] + C7 fjf |0 fc - ^| 2 dW"" 1 ) 

CO. 

Similarly we obtain 



I m Jm 
— > as fc — > oo. 



G(x,$ fc (y))^(x)^ fc (2/)J M $ fc (2/)^"- 1 (x)dH"- 1 (2/)^ / / G(x, yMx)ip(y) dU n ~ Y (x)dH n - 1 (y). 

M J M JmJm 

These two estimates yield (|5.33[) and the lemma is proved. 

□ 



We are now ready to prove Proposition 15.21 

Proof of the Proposition 15.21 Without loss of generality we may assume that F is a G 2 regular and 
Mf = dF n il is a G 2 -manifold with boundary and it meets dfl orthogonally Moreover we assume that 
\\F, E\\\y2, P < S where 6 > is chosen later. 

Suppose that X is the vector field and $ is the associated flow given by Lemma 15.31 For every 



t £ [0, 1] we denote E t = &(E, t), M t = dE t fl fl and v t as the unit normal of E t . Since E is a critical set 
we have 

Moreover Lemma ISTTI yields ^ J(£ l t )| t _ Q > ^-(MI^m) where (p = (X, v) for some constant c > 0. The 
idea of the proof is to show by a continuity argument that when 5 > is chosen small enough we have 

(5-34) ^J(E t ) > ^\\<Pt\\ 2 HHMt ), 

for every t <G [0, 1] where ip t = {X, u t ). 

The result follows from the estimate (|5.34j) . Indeed when 8 > is small we have by (|5.28|) and (|5.30l) 

that 

\ y > t \*dU n - 1 =( J M § t \(X,(y t o$> t ))\*dH n - l >c [ \{X, v )\ 2 dU n -\ 
M t Jm Jm 
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Therefore the estimates (|5.34|) and (|5.26j) from Lemma T5. 41 yield 

J(F) - J(E) = J (1 - t )^J(E t )dt > c\Ml HM) >c\FAE\ 2 , 

and the claim follows. Hence, we need to proof (|5,34[) . 

The formula for ^ J(E t ) in Proposition 14. 1 1 can be written as 

d2j }^ f ) = d 2 J{E t )[ Vt ] - [ (H 9Et + A 1VEt ) div Tt (X Tt (X, v t )) dK 1 - 1 

+ [ (H dEt + 4 7 «E t ) div(X)(X, v t ) dn n -\ 

where d 2 J(E t )[] is the quadratic form defined in (|2.6p and (p t = (X,v t ) as before. Notice that since the 
flow preserves the volume we have 



f cptdH 11 - 1 = f {X,v t )dU n ~ l =Q. 

JMt JMt 



When 5 > is small enough Lemma [531 yields 

(5-35) d 2 J(E t )& t ] > |||^||| 1(Mt) . 

We proceed by showing that when S > is small enough we have 
(5-36) S t := | / (H dEt + 4 1VEt )div Tt (X Tt ((X,v t ))) dK"- 1 ] < C -^\Wt\\ 2 HHMt y 

Indeed, since E solves (|2.4j) we deduce (Hg Et + 47^) o — > A in L P (M). Therefore when S > is 

small we have for small e > 



St = I / (H dFt + 4 7 «F t - A) div rt (X Tt (X, v t )) dU n - 

JM t 



< ( / (H m +4 7 v Ft - Xfdn n - L / (div Tt (X Tt (X,v t )))— dU r 
hit J \JM t 



< e (\\D Tt X Tt \\ LHMt) \\(X,v t )\\ j. + \\D Tt (X,v t )\\ L 2 {Mt) \\X 

\ Lp- 2 (Alt) 

<Ce\\D Tt X Tt \\ m{Mt) \\(X,v t )\\ m(Mt} 



Tt I I 2 P 

Lp- 2 (M t ) 



where the last inequality follows from Sobolev inequality and the fact that p > n. The estimate (|5.36[) 
then follows from (|5.27[) . 

Similarly we prove that when S > is small we obtain 

(5.37) R t := I / (H dFt +4 1 v Ft )div(X)(X,v t )dH n - 1 \ < ^||^|| 2 ffl(Mt) . 

The volume constraint \F t \ = \E\ yields 

dw(x)(x,v t )dn n - 1 =0 



'M t 

as it was discussed earlier after Proposition 14.11 As before using the fact that E solves (|2.4[) we obtain 
by Holder and Sobolev inequalities 



Rt = I / (H dFt + 4~/v Ft - A) div(X)(X, v t ) dK 1 - 1 1 

JM t 

<( f ( H dFt + ^v Ft - A) p dH' 1 - 1 ) P ( f (div(X)(X, v t ))^ dW 

\JM t J \JM t 

<Ce\\dw(X)\\ LHMt) \\(X,v t )\\ Hl{Mt) , 
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for any e > when S > is small enough. The estimate (|5.37|) follows from (|5.27l) . 

The key estimate (|5.34|) now follows from (|5.35|) . (|5.36j) and (|5.37|) and the proposition is proved. 

□ 



6. Proof of the main theorem 



In this section we prove Theorem 12.31 As it was mentioned in the introduction we will effectively 
use the regularity of u-minimizers to rule out those competing sets which are not regular. This idea goes 
back to the work by Cicalese and Leonardi [8] and by Fusco and Morini [9]. 

We begin with a simple lemma. 



Lemma 6.1. Suppose that E C ft is C 2 -regular and satisfies the orthogonality condition (|2.5[) . There is 
a constant C depending only on E such that for every F <S BV(fl) we have 

J{E) < J{F) + C\FAE\ . 

Proof. By the Lipschitz continuity of the non-local part (|3.2j) we have 

J{E) < J{F) + P Q (E) - P Q (F) + C \FAE\ . 

On the other hand by the assumptions on E we may construct a C 1 vector field X on ft such that X = ve 
on M = dEnfl, (X, v n ) = on dfl and \\X\\ L oo < 1. Therefore 

Pn(E)-Pn(F)< [ (X,u) dU n -\x) - f (A» dU n -\x) 
Jd(Enn) Jd(Fnn) 

= / div X dx — / div X dx 

JEnn J Fnn 

< ||divX|| L oo|FA£;| . 

□ 



We remark that the above proof also yields that E is an w-minimizer in the sense of Definition 13.41 
The next lemma states that the convergence of mean curvatures implies W 2 p convergence for any 
p > 1. It is very similar to Lemma 7.2 in [T]. 

Lemma 6.2. Letp > 1. Suppose that E is as in Theorem \2.3\ and there are diffeomorphisms $>k ■ E — > F^ 
such that \\$k — Id\\ci,<* — > and 

H dFk ($k(-)) -> H dE (-) in If. 
Then there are diffeomorphisms $fc : E — > F^ such that ||$fc — /dH^.p — > 0. 

Sketch of the proof. Since F^ are C 1,a manifolds with boundary they can be presented locally as a graph 
of C > a functions. Away from the relative boundary dFk n il the convergence follows as in (pQ, Lemma 
7.2). On the boundary we follow a standard argument and flatten the boundary dil locally to half space 
and then use similar elliptic regularity estimate as in the interior case. □ 

We are now ready to prove the main result of the paper. As in [TJ E] the proof is divided into two 
steps where we first prove that E is a local minimizer with respect to Hausdorff distance and then sharpen 
the result to obtain the local minimality with respect to L 1 distance with the quantitative estimate. 
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Proof of Theorem 1 2. 3[ Step 1: First we prove that E is a local minimizcr in L°°-topology, i.e., there 
is 5 > such that for every F 6 BV(Q) with \F\ = \E\ and Fc{ieO dist(x, E) < 5} it holds 

J(F) > J(E). 

For simplicity we denote ls{E) := {.t 6 £7 | dist(a;, E) < 6} 

We argue by contradiction and assume there are Ek C f2 such that \Ek\ = \E\ , Ek C Ii/^{E) and 

J(E k ) < J(E). 

We will use the signed distance function dn,£;( - ) of E on fi constructed in the first step of Lemma [5.31 
Denote the sets {x <E il \ dn^ix) < 1/fc} U E by U k . The reader needs only to know that these sets are 
C 3 regular and satisfy 

(i) T 1 / k ,(E) C U k , for some k' 

(ii) ||{7fc,S||c2 —> in a sense of definition ()5.1|) and 

(iii) Uk satisfies the orthogonality condition (|2.5j) . 

As in [I] we replace the contradicting sequence (Ek) by (Fk) which minimizes 

(6.1) J(F) with constraint \F\ = \E\ and F C f/ fc . 

By the contradiction assumption J(-Ffc) < J(E) for every A: and obviously Fk ~ > S in Hausdorff topology. 

Let us show that -Ffc are w-minimizers with uniform constants A and r. To that aim let G C be a 
regular set such that GAFk C B r (xo). Let us divide G into two parts 

G(lUk and GUU k . 

Since minimizes (|1.3[) in C/ fc with a volume constraint we obtain as in the proof of Proposition l3.4l that 

(6.2) PB r ( Xo )(F k ) < P Br{xo] (GnU k ) + C\(GAF k ) n t/ fc | 

for some C. Moreover, since Uk are uniformly C 3 regular and satisfy the orthogonality condition they 
are w-minimizers, as we remarked after the proof of Lemma l6.1[ and therefore 

(6.3) P BrM (Uk) < P Br{xo) (GUU k ) + C\G\U k \, 
for some C > 0. Since 

p Br(xo) (Gnu k ) + p Br(xo) (GnUk) < P Br{xa) (G) + P Br{xo) (u k ) 

the estimates (|6.2[) and (|6.3p yield 

^(, )(^) < P Br{xo) (G) + A\F k AG\ 

for some large Aand the cj-minimality is proved. By the w-minimizing property of Fk and by Theorem 
13.51 we conclude that Fk — > E in C 1 '" and that Fk satisfy the orthogonality condition (|2.5[) . 

Denote the relative boundaries as usual by M = dE n O and Mfc = dF k fl fi. The Euler-Lagrange 
equation for F/. reads as 

J #SF fc + 47ii Ffc = X k on M fc n E/ fe , 

(6.4) i 

[ i/aF fc + 4:-fv E + lk = A on M*. \ C/ fe , 

where 7^ are some remainder terms with 7^ — > as k ~ > 00. Since Ft are C 1 regular Proposition 13.21 
yields ||HaF fc ||L=° < A. Moreover the equation (|1.4[) implies 

(6.5) -> in C 1 ^). 
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Next we show that 
(6.6) H dFk ($ k (-)) ^ H aE (-) oni"(M), 

for any p > 1. Indeed, consider a vector field X S C 1 (r2,R n ) as in the proof of Lemma UTTl i.e., X = vm 
on M and (X, i/q) = on dCl. We multiply the equation (|6.4[) by (X, ^j\/ fc ), integrate over Mfc and then 
integrate by parts and use the C 1, "-convergence to deduce 



/ {\ k -A 1 v Fk ){x, VMk )<m n - 1 + I {\-A 1VE - lk ){x, VM )dn 

JM k nu k JM k \u k 



n-1 



[ H dFk (X,v Fk ) dH"- 1 = f div Mk XdH n ~ 
J M k nn J M k nn 



-> / div M X dH 
'Mnn 



n-1 



{\-^v E )dn n -\ 

Mnn 

where the last equality follows from the Euler-Lagrange equation (|2.4I) for E and from X = v E on dEDtl. 
The C 1,a -convergence of and (|6.5p then imply that either A^ — ^ A, or H n ~ 1 (M k \ U k ) 0. In cither 
case we obtain (|6.6[) since Hg Fk are uniformly bounded in L°° . 

From (|6.6[) and Lemma [6~2l we deduce ||Ffc,£J|| W 2, P — ► 0. Since F k satisfy the orthogonality condition 
we may use Proposition 15 . 21 to conclude 

J{F k ) > J(E) 

when k is large. This contradicts with the minimality of F k (|6 . 1 1) since by the contradiction assumption 
we have J(F k ) < J{E). 

Step 2: We are ready to prove Theorem 12.31 

As before we argue by contradiction and assume that there are E k C SI such that \E k \ = \E\, 
\E k AE\ -> and 

J(E k ) < J(E) + ^ \E k AE\ 2 , 

for ci > from Proposition 15.21 Denote e k := \E k AE\. As in [T] we replace the contradicting sequence 
(E k ) by (F k ) which minimizes 



(6.7) min{ J(F) + C ^ (\FAE\ - e k ) 2 + e k :FcS! with |F| - \E\} , 

for some constant C which will be chosen later. 

By compactness we may assume that, up to a subsequence, F k — > Fq in L 1 and that Fq solves 

mm{ J(F) + C\FAE\ : F C fi with \F\ = \E\] . 

By choosing C large enough it follows from Lemma I57T1 that F = E. In particular F k —> E in L 1 . 

As in step 1 we observe that every F k is an w-minimizer with uniform constants A and ro- In fact, 
since there are no obstacle in (|6.ip this observation follows exactly as in the proof of the Proposition ^. 41 
Therefore Theorem 13.51 implies that F k — > E in C 1 '" and that F k are C 1, "-manifolds with boundary for 
k large and satisfies the orthogonality condition (|2.5ll . 

The Euler-Lagrange equation for F k reads as 

f H dFk + Ajv Fk + C jMMk^iL = A fe on M k \ M, 

(6.8) { y/(\FAE\-e k ) 2 +e k 

{ H dFk + 4~/v E = A on M k n M. 

Since the mean curvature Hd Fk is bounded we conclude that M k is W^ ,p -regular. 
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Using the minimality of F k and the contradiction assumption and step 1 we obtain 
J(F k ) + C^(\F k AE\-e k ) 2 +e k < J{E k ) + CJT k 



<j( E ) + c fel + cvr k 

< J(F k ) + ^e 2 k + C^T k , 

where the last inequality follows from the L°°-local minimality proved in step 1. The previous inequality 
yields 

(6.9) Urn mm-**) = i 

fc-s-oo Eh. 



In particular 



hm , (mSLM = 0. 



y/(\FAE\~e k ) 2 +e k 

Arguing exactly as in step 1 we conclude from the equation (|6.8[) and from (|6.Q[) that | \F k , E\ \ w^.p — > 
for any p > 1. Hence we may use Proposition 15 . 21 to conclude 

J(F fe ) > J(E)+ Cl \F k AE\ 2 . 

However the minimality of F k , the contradiction assumption and (|6.9[) yield 

J(F k ) < J(E k ) < J(E) + C ±e\ < J(E) + j\F k AE\ 2 

which is a contradiction. □ 

Proof of the Corollary \2-4\ We obtain Corollary 12.41 immediately from Theorem 12.31 arguing as in 
Theorem 6.3 in pQ. □ 
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